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A method of estimating the volume flow in a 
vessel utilizing ultrasound techniques including 
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ultrasonic transducer array so as to produce fluid 
velocity information for the vessel, the plane 
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having a direction which is offset from an axis of 
the vessel such that the beam plane forms an 
elliptical section within the vessel, and utilizing 
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mean velocity for the fluid flow through the 
elliptical section. The improvement of the volume 
flow estimation is achieved by utilizing average 
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and mean velocity calculation based on a 
weighted summation of the velocity values a 
multiple measurement points with weighting 
factors determined by the Doppler signal power 
at the measurement points. 
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2.**** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 



CLAIMS 



[Claim(s)] 

[Claim 1] It is the approach of evaluating the volume flow in tubing using an ultrasonic technique. The 
ultrasonic transducer array which generates a beam so that fluid rate information may be made from tubing 
is used. Are the step which picturizes tubing and the field formed of two or more beams generated by the 
ultrasonic transducer The step to which it has the direction [ shaft / of tubing ] shifted and an ultrasonic 
transducer array picturizes the ellipse cross section in tubing, step which calculates the mean velocity about 
the fluid which passes along an ellipse cross section using the aforementioned rate information about two or 
more points in tubing and an ellipse cross section flowing from — the approach of becoming. 
[Claim 2] Furthermore, the approach characterized by consisting of a step which doubles an ellipse cross 
section with the aforementioned fluid rate information, and a step which determines the ellipse cross-section 
parameter for tubing from the aforementioned ellipse cross section with which it doubled. 
[Claim 3] It is the approach characterized by the aforementioned tubing being a man or one part of the 
inside of the body of an animal in the approach indicated by claims 1 or 2. 

[Claim 4] It is the approach characterized by the aforementioned approach being performed by the 
reprogramming of the software of a standard ultrasonic device by the aforementioned ultrasonic device in 
the approach indicated by either of claims 1-3. 

[Claim 5] The step which calculates the aforementioned mean velocity in the approach indicated by either of 
claims 1-5 is an approach characterized by relating the weighting **** average with the magnitude of the 
input signal in the point corresponding to the aforementioned rate information by this power including using 
the weighting **** average by the power of the aforementioned rate information. 

[Claim 6] The approach characterized by processing the aforementioned rate information beforehand in the 

approach indicated by either of claims 1-5 so that the rate data value in the area of the echo of the effective 

(significant) gray scale from the aforementioned ultrasonic transducer may be made small. 

[Claim 7] The approach characterized by throwing away weighting of power smaller than the magnitude 

decided beforehand in the approach indicated by claim 5. 

[Claim 8] Equipment which enforces one approach of claims 1-7. 



[Translation done.] 
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DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[Field of the Invention] 

This invention relates to utilization of the field of ultrasonic imaging technique, and the ultrasonic data 

which create the information on a volume flow especially. 

[0002] 

[Description of the Prior Art] 

In recent years, an ultrasonic image pick-up is becoming still more important in the field of a medical- 
application image pick-up of a man or the inside part of the body of an animal especially. The principle of 
an ultrasonic image pick-up is known well. However, the activity of the phase data from an ultrasonic signal 
flows in relation to an ultrasonic scan, and is creating information recently. Many manufacturers offer the 
"color Doppler image pick-up" which offers information in the form of a color guide (indicator) which 
shows the direction and rate of the rate of flow of blood in tubing by the end of today. 
[0003] 

[Problem(s) to be Solved by the Invention] 

Ideally, to use the hardware of standard marketing only for software, changing is desired for noninvasive 
real-time supersonic-wave measurement of the volume rate of flow of the blood in tubing in the certainty of 
high degree of accuracy and high extent. 
[0004] 

The object of this invention is offering exact measurement of the flow volume in tubing using ultrasonic 

technique. 

[0005] 

[Means for Solving the Problem] 

From the 1 st standpoint concerning this invention, the approach of evaluating the flow of the volume in 
tubing using an ultrasonic technique The ultrasonic transducer array which generates a beam is used so that 
fluid rate information may be made from tubing. Are the step which picturizes tubing and the field formed 
of two or more beams generated by the ultrasonic transducer It has the direction [ shaft / of tubing ] shifted 
and consists of a step which calculates the mean velocity about the step to which an ultrasonic transducer 
array picturizes the ellipse cross section in tubing, and the fluid which passes along an ellipse cross section 
using rate information about two or more points in tubing and an ellipse cross section flowing. 
Preferably, mean velocity is calculated using the weighting **** average by the power of the 
aforementioned rate information, and the weighting **** average is related with the magnitude of the input 
signal in the point corresponding to the aforementioned rate information by this power. 
It consists of a step which doubles an ellipse cross section with the desirable fluid rate information of further 
the above [ the aforementioned approach ], and a step which determines the ellipse cross-section parameter 
for tubing from the aforementioned ellipse cross-section parameter with which it doubled. 
The aforementioned approach suits ideal for the activity about a man or one part of the inside of the body of 
an animal, and can be performed by the ultrasonic device by the reprogramming of the software of a 
standard ultrasonic device. 
[0006] 

[Embodiment of the Invention] 

In the gestalt of desirable implementation of invention, in order that the slanting scanning may extract the 
detail of the ellipse cross-section structure of tubing, and in order to use an elliptical ring structure, it 
combines with the fluid migration information in the cross-section structure, and is used so that a volume 
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flow may be drawn. 
[0007] 

The usual sonication technique is first shown in drawing 1 in graph. In this usual arrangement, a transducer 
1 generates the beam 2 of a plane configuration substantially [ the thickness decided beforehand ]. This 
beam 2 is processed so that it may be projected through tubing 3, that echo may return, and it may be 
received by the transducer 1 , then the structure of tubing may be determined. In the color Doppler 
picturizing method, generating of a series of supersonic waves is performed continuously quickly 
periodically, and the echo of return is recorded. The phase contrast between multiplex return is determined, 
and it is used in order to calculate the rate-of-flow information in tubing 3 by the known approach. 
[0008] 

A slanting scan is used in the gestalt of desirable operation. Here, bearing of a beam 2 is not meeting tubing 
1 . As especially shown in drawing 2 and drawing 3 , the slanting scan beam revolution of two forms is 
given. As shown [ 1 st ] in drawing 2 , the 1 st cross-section volume 4 is produced using the revolution alpha 
of a surrounding probe (namely, surroundings of the longitudinal direction shaft of a transducer probe) of 
the shaft 7 vertical to the shaft of the longitudinal direction of tubing. As shown [ 2nd ] in drawing 3 , the 
surrounding revolution beta of the shaft 8 vertical to a scanning field is given. In this way, the final cross 
section 4 of the obtained tubing 3 is the configuration of an ellipse cross section, and this is volume which is 
determined with the thickness property of a transducer beam and which was decided beforehand. 
[0009] 

in the gestalt of desirable operation, a series of frames picturize using the slanting scan which is shown in 
drawing 2 and drawing 3 and which makes a certain include angle — having — next, the target tubing — ** — 
in order to determine the actual volume rate of flow to break, it is processed so that it may explain later. 
[0010] 

In drawing 8 , two or more steps 1 0 which can be set to the approach of the gestalt desirable operation are 
shown. By this approach, an ultrasonic facility of the present age of the criterion which picturizes a series of 
frames using the slanting scanning is first used at step 11. 
[0011] 

Corresponding rate information is acquired from a series of picturized frames using a well-known color 
Doppler technique (step 12). The method of obtaining color flow data is completely standard. 
[0012] 

Next, the image of color flow is examined so that an ellipse cross section may be identified, and the shape 
parameter (alpha and beta) of a scan is obtained after this (step 13). The ratio of the minor axis and major 
axis of an ellipse is used in order to determine an include angle alpha, and on the other hand, in order that 
the direction of the major axis over a horizontal line and a color Doppler setting-out angle may determine an 
include angle beta in a scan, it can be used. 
[0013] 

If the whole elliptical ring structure is determined next, it will be used in order to calculate the cross section 
of tubing. Next, the estimate of the volume rate of flow can be calculated by having the information of a 
configuration and mean velocity (step 14). 
[0014] 

A more detailed understanding of this processing can be acquired by examining the step used in order to 

simulate a scan and signal processing mathematically. 

[0015] 

In order to determine and identify an ellipse parameter, and in order to calculate a volume style, the image 
pick-up side received by the transducer component is scanned continuously, and the tubing coordinate 
which corresponds about each scanning coordinate is acquired for processing which determines whether a 
sample coordinate (this expresses the sample voxel) is in the lumen of tubing. 
[0016] 

Therefore, it is required to change between scanning field coordinate space and tubing coordinate space. 
This transformation method is explained with reference to drawing 6 from drawing 4 . In drawing 4 , x of 
tubing, y, and a z-coordinate system are shown. Drawing 4 shows the transducer 1 located in a line in the z- 
y side again. Although a tubing coordinate remains in x, y, and z criteria frame in this include-angle slanting 
method, a scanning field coordinate is changed into x" and y"z n criteria frame using two revolutions. 
[0017] 

As shown in drawing 5 , the 1st revolution is a revolution of the surrounding include angle alpha of the z- 
axis shown in drawing 5 by the transducer 1 . If the notation of a matrix is used, this conversion can be 

http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_webcgi_ejje 1/23/2006 



Page 3 of 6 



written as follows. 
[Equation 1] 
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[0018] 

Direction cosine is as follows from drawing 5 using this notation. 
[Equation 2] 
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Namely, [Equation 3] 
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[0019] 

the 2nd conversion shown in drawing 6 ~ setting — the revolution of a transducer — x" from x'-y'-z' frame of 
reference ~ it changes to -y"-z" frame of reference. This conversion is expressed as follows. 
[Equation 4] 

(cosZ(x',x") cosZ(y',x") cos Z(z' ,x")Y x'} 



r 



K Z J 



cosZ(jc',y) cos Z(y',y") cos Z(z',y") y' 
^cosZ(x',z") cosZ(/,z") cosZ^z")^', 



(4) 



[0020] 

It is as follows using the notation of drawing 6 about direction cosine. 
[Equation 5] 

r cos(0) cos(*r/2) cos(tt/2) 



y 



cos{n 12) cos( J3) cos(tt 1 2 + p) 
cos(;r 12) cos(tt / 2 - J3) cos(/?) 



y 



(5) 



Namely, [Equation 6] 

0 0 

0 cos(a) - sin(>#) 

^0 +sin(yff) cos(^) J 



y 



K z j 



^x^ 



k z s 



(6) 



[0021] 

Two conversion is combined and one conversion to x"-y"-z" frame of reference from the x-y-z frame of 
reference is obtained from a formula (3) and a formula (6). 
[Equation 7] 

0 0 Y cos («) -sin(a) OYx^ 

sin(a) cos(a) 0 y (7) 
0 0 1 



X Z j 



0 cos(a) - sin(/?) 
0 sin(yfl) cos(>9) 



y 

A Z J 



http://www4.ipdl .ncipi . go.jp/cgi-bin/tran_web_cgi_ejj e 



1/23/2006 



JP,2002-505901,A [DETAILED DESCRIPTION] 



Page 4 of 6 



From now on, the following formula will be obtained. 
[Equation 8] 

f cos(a) -sin(a) 0 

cos(>#) sin(a) cos(/?) cos(a) - sxn(/3) 
k sin(/?) sin(a) sm(j3) cos(a) cos(/?) 



\ Z J 



J\ Z J 



(8) 



f x\ ( cos(a) 



Or it is [Equation 9] by matrix conversion. 

cos(y?) sin(a) sin(/?)sin(a) 
- sin(a) cos(/?) cos(a) sin(/?) cos(a) 
0 -sin(>9) cos(>£) 



\ Z J 



(9) 



The matrix of a formula (8) and a formula (9) is mutual rewriting at an interesting thing. 
[0022] 

In this simulation "idealized", the increment of the location (0, y", z") of the sample volume is carried out 
over the breadth of a scanning field, and the coordinate (x y, z) of corresponding tubing is calculated using a 
formula (9). Supposing the sample volume is in the interior of the lumen of tubing (namely, [Equation 10]) 



x 2 +y 2 



< Radius 



The value of the power which **** and corresponds is assigned to the value of 1 .0, and a rate is assigned to 

the value based on a distance radial [ from a tube axis ]. 

[0023] 

Though regrettable, it is disregarded that the beam of a transducer has the width of face of finite, and above- 
mentioned analysis has the cross section of Gaussian distribution on the strength mostly. Therefore, in the 
power in the location of each sample volume, i.e., the voxel, and count of a rate, it is dramatically desirable 
to include an ultrasonic parameter. 
[0024] 

Decision of a configuration (ellipse doubling) (step 13 of drawing 8 ) 

The 1st step of this approach is the decision of a configuration, and this requires doubling to the power data 
or rate data of an ellipse. The data used in order to double an ellipse include the combination of power data, 
rate data, and gray scale information. In one example, power data are used in order to double an ellipse, they 
hold down a power value smaller than a threshold, and make small the field which has an effective 
(significant) gray-scale value. 
[0025] 

The algorithm used for ellipse doubling is based on the formula of moment of inertia. About rate data, only 
the value beyond a predetermined rate threshold is used and each value has equal weight in count of moment 
of inertia. On the other hand, when setting power data by the ellipse, the value beyond a predetermined 
power threshold is used as weight in count of moment of inertia. In this way, more fields to count of higher 
power than the field of lower power contribute. 
[0026] 

If it doubles with an ellipse, an include angle alpha and an include angle beta are calculable as follows. 
Alpha=sin -1 (b/a) 

It is a here. = The die length of the major axis of an ellipse b = The die length of the minor axis of an ellipse 
= Diameter of tubing beta = Include angle between the major axes of an ellipse and horizontal lines in a scan 
[0027] 

From alpha and beta, "Doppler include-angle" amendment can calculate as follows. 
[Equation 11] 

C °S(& Doppler ) = COS(a)COS(£> - /?) 

Although psi is a color Doppler setting-out angle and is 70 degrees, 90 degrees, or 1 10 degrees typically, it 
is not limited here by it. The relation between a, b, beta, and psi is calculable as follows using an include 
angle alpha. 
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[0028] 

(It is vertical to the shaft of a blood vessel) The true area Aactual of a cross section is calculable as follows 
from the measuring-plane product Ameasured using an include angle alpha. 
[Equation 12] 

Actual = Measured sin(a) (14) 

[0029] 

Count of mean velocity (step 14 of drawing 8 ) 

As stated above, when a shape parameter is measured from an ellipse cross section, the following step in a 
desirable operation gestalt is calculating mean velocity covering the volume divided by the ellipse. 
[0030] 

Preferably, eclipse ******** with weight is used by the power described in a following formula. 
[Equation 13] 

v = — (15) 

(y" 1 z n )wthin ellipse 



P expresses the received power value here and nu expresses the received rate value, moreover (y", z") It is in 

an ellipse. 

[0031] 

This way gives small weight to the voxel of lower power, for example, the voxel which accepts it selectively 
and is in tubing. Therefore, ********** approach turns caution to the "subfield product" effectiveness 
relevant to the voxel dimension of finite by this power. 
[0032] 

Of course, angle correction is an average reading per second [several 14]. 



* true 



since — true mean velocity [several 15] 



y 

measured 



It is required for a ****** sake. Such amendment is applicable as follows. 
[Equation 16] 

Vtrue = x (1 6) 

c °s(0 Wr ) 

[0033] 

Estimate of the volume rate of flow (step 14 of drawing 8 ) If shape parameters alpha and beta, the measured 
area, and mean velocity are given, the volume rate of flow Q is calculable as follows using a formula (13), 
(14), and (16). 
[Equation 17] 

Q=v-A= V ™™ d 

COStf^) (i7) 

_ U measured ^measured S ^ n (^) _ V measured ^measured ) 

cos(a ) cos((p - J3) cos(#> - p) 

[0034] 

Since it stated above, it turns out that it flows with a rate and the more exact result of the volume is obtained 

using a power weighting system. 

[0035] 

Of course, many deformation is possible. Such one deformation is using "expected power weighting" of 
modeled power weighting about a certain ellipse cross-section structure. This makes light the problem 
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relevant to power data with many noises which may arise. 
[0036] 

Another deformation is being able to use it in order to process flow information, in order that location data 
may photo an organization echo etc. This deformation includes the activity of this information for holding 
down the rate data value in the area of an effective (significant) organization echo by the trial of the 
traditional gray-scale data usually given to an ultrasonic device, and fake data. 
[0037] 

Other many deformation is possible, without separating from the summary of this invention which indicated 
widely this invention shown in the above-mentioned specific operation gestalt so that I may be understood 
by this contractor. Therefore, it should be understood that this operation gestalt is a thing for explanation of 
invention, and is not restrictive. 
[Brief Description of the Drawings] 

[Drawing 1 ] Drawing showing a standard longitudinal direction scan [Drawing 2] Drawing showing the 2nd 
revolution of a scanning field relative to tubing [Drawing 3] Drawing showing the 2nd revolution of a 
scanning field relative to tubing measured [Drawing 4] Drawing showing the process of coordinate 
transformation of being used in a desirable operation gestalt [Drawing 5] Drawing showing the process of 
coordinate transformation of being used in a desirable operation gestalt [Drawing 6] Drawing showing the 
process of coordinate transformation of being used in a desirable operation gestalt [Drawing 7] Drawing 
showing the relation between the variables used in a desirable operation gestalt [Drawing 8] Drawing 
showing the step in a desirable operation gestalt [explanation of a sign] 
1 Transducer. 



[Translation done.] 
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* NOTICES * 

JPO and NCIPI are not responsible for any 
damages caused by the use of this translation. 

1 .This document has been translated by computer. So the translation may not reflect the original precisely. 
2.**** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 



DRAWINGS 



[Drawing 1] 




[Drawing 3] 
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[Drawing 4] 
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z" -y" BHzfclt* h7>Xfa- 



ylft;*fi<Z>lfimtt 



[Drawing 71 
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[Drawing 8] 
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[Translation done.] 
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